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Abstract — We calculate the electronic and thermal transport 
properties of devices based on finite graphene antidot lattices 
(GALs) connected to perfect graphene leads. We use an atomistic 
approach based on the 7r-tight-binding model, the Brenner poten- 
tial, and employing recursive Green's functions. We consider the 
effect of random disorder on the electronic and thermal transport 
properties, and examine the potential gain of thermoelectric 
merit by tailoring of the disorder. We propose several routes to 
optimize the transport properties of the GAL systems. Finally, 
we illustrate how quantum thermal transport can be addressed 
by molecular dynamics simulations, and compare to the Green's 
function results for the GAL systems in the ballistic limit. 

I. Introduction 

Thermal management in nanoelectronics is a topic of great 
importance! 1 1, and it is relevant to consider materials for 
thermo-electric cooling at the nanoscale. Although graphene 
is an excellent heat conductor, its thermoelectric properties 
seem to be highly tunable in nano structured graphene, such 
as graphene antidot lattices (GALs) (2). GALs have been 
proposed as a flexible platform for creating a semiconducting 
material! 3], E), Q, O and the regular nanoperf oration of a 
graphene layer may be a solution for making graphene a ver- 
satile material for electronics as well as thermal management. 

Several groups have developed fabrication techniques for 
creating this nano structured system, e.g. by electron beam 
lithography 171. 151 or block copolymer lithography 0, ifTOl . 
The latter method results in a highly regular hexagonal lattice 
of holes. In another experiment hexagonal holes were obtained 
with very pure edge chirality by anisotropic etching (TH. 
Various types of disorder play an important role for trans- 
port properties. Edge disorder in nanoribbons has previously 
been theoretically predicted to increase the thermoelectric 
efficiency C21 and widen the band gap distribution of GALs 
with periodic disorder fT3l . Recent experiments have managed 
to control the concentration of the 13 C isotope in both amount 
and space |14'| and it has also been possible to measure 
the thermal conductivity by Raman spectroscopy fT5l Several 
theoretical predictions have been made on the isotope effect 
in graphene (16]|, ifTTl . nanoribbons fT8lL fT9lL l20l and carbon 
nanotubes l2n . 

Here we use atomistic modelling to address how structural 
and isotope disorder can reduce the thermal conductance of 
GAL devices. 



II. The role of disorder in GALs 

The electronic and phonon transport properties are cal- 
culated from the energy-dependent transmission functions, 
T e and T p h, using a Landauer-type formula. The Landauer 
formula for the electrons reads 
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where uf(E, Vl/r) * s me Fermi-Dirac distribution at the 
chemical potential of the left/right lead. Within the linear- 
response limit we can consider variations with changes in the 
chemical potential, /i, e.g. by doping or gating of the graphene 
system, by evaluated the following moments of the electronic 
transmission, 
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They relate the electronic current and the electron heat current, 
in the linear response regime. Thus, we have (22) the electrical 
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takes an analogous form and within linear response the thermal 
conductance is given by 



e 2 Lo, the electron thermal con- 
/T, and the Seebeck coefficient 
. For phonons the Landauer formula 



f 

JO 



du 



\%h{u)- 



e * B ± 



(3) 



From these we can evaluate the thermoelectric figure of merit, 
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We consider only disorder in the transport direction in GAL 
structures with periodicity in the direction perpendicular to 
transport, see Fig.[T} and employ averaging of the transmission 
functions over a dense k± -point sampling. Using the Brenner 
potential O we perform structural relaxation, obtain the dy- 
namical matrix, and then calculate the electronic Hamiltonian 
using the 7r-tight-binding model for the device region of a 
given length between pristine graphene leads. Based on this we 
calculate the /c-averaged T e and T p h by an atomistic Green's 
function method, which is efficient for long systems with 
different kinds of disorder l24l . l25l . Note that we assume the 
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Fig. 1. System setup and the computational rectangular unit cell (green 
rectangle). Two graphene leads are connected by the finite GAL of length 
2 (M = 2) corresponding to 4 holes in the direction of transport. The 
depicted system is based on a {L = 10, S = 3arm} GAL, where we use 
the nomenclature from |2|. The red atoms are randomly removed during the 
generation of edge noise. The unit cell is 10 four- atom armchair units wide. 
We apply periodic boundary conditions in the transverse direction. 
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device to be free-standing, and neglect all effects of electron- 
phonon interaction such as phonon-drag. Later, we comment 
on how anharmonic effects may be included using molecular 
dynamics. In the remaining of this section we apply this 
method to three types of disorder; edge disorder, position 
disorder, and the isotope effects. 

A. Lattice disorder 

In Fig. [2] we show the outcome of transport simulations 
for various GALs with random edge disorder. The system is 
constructed from a sequence of GAL unit cells (with two 
holes) in which a chosen number of atoms are randomly 
removed at the edge one by one. We start from a very small 
{10,2.3arm}-GAL, and increase the hole size in this random 
fashion until it is comparable to the {10,6arm}, which have 
a max(ZT) « 0.25 at 300K[2J. Here, max(ZT) refers 
to the maximal occurring ZT within a gate bias range of 
[— 1.5; 1.5]eV. We do not allow for C-atoms with a single 
nearest neighbor, since these are energetically less favorable. 
Regarding edge disorder the decrease of thermal conductance 
outweights the decrease in the power factor (S 2 G e ), as can 
be seen in the increase in ZT with the number of removed 
atoms in Fig. [2] However, the effect for disordered antidots is 
much smaller as compared to the effect for hole-size for edges 
with pure chiralityEl. Therefore, one should preferable avoid 
the edge disorder during fabrication if one seeks to optimize 
thermoelectric applications. We note that the kinetics of the 
system or electronic current|26] may stabilize pure edges to 
some extent. 

Next, we examine disorder in the hole lattice positions by 
randomly shifting a sequence of 2-hole unit-cells along the 
transverse direction by dN (armchair unit periodicity). The 
resulting figure of merit is given in Fig. [3] The full cell is 
L = 10 armchair units wide, and we examine dN = 1, 2, ..5 
corresponding to increasing lateral disorder. The electronic 
conductance is almost unchanged at a low degree of position 
disorder (dN = 1), and the power factor can even initially in- 



Fig. 2. Edge disorder in a GAL system with 5 unit cells along the device, 
corresponding to a length of 380A. The system is based on a {10,2.3arm} 
GAL (with a width of 43A) and we increase the number of atoms removed 
at random positions at the edge (the red atoms in the insert, showing one unit 
cell at N rern /Ntot = 0.19). The ratio of removed atoms to unit cell size 
is changed through the interval previously considered for perfect antidots | '2' |. 
The resulting edges are highly disordered and can be taken as a worst case 
scenario (e.g. by vapor deposition of carbon atoms on a surface). 
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Fig. 3. Position disorder is here studied by increasing an interface mismatch 
between two GAL regions. Increasing the position disorder tends to reduce 
both the thermal conductance and also the power factor except at a very low 
degree of disorder. The system has a length of 380A(M = 5) and is based 
on a {10, 6arm} GAL. 

crease at a low degree of position disorder, before it eventually 
decreases. Meanwhile, the thermal conductance decreases with 
increasing dN. Therefore, especially a weak degree of position 
disorder increases ZT to above 0.5. For perfect GALs with 
a high antidot density the electronic ZT (setting K p h = 0) is 
not much larger at the chemical potential of the largest peak. 
For low /jl the n p h dominates, and therefore it is the first ZT 
peak which is increased by the disorder. On the other hand, the 
highest peak at higher \i mainly follows the change in power 
factor. 

B. Isotope disorder 

Now we address the effect of isotopes positioned at random 
sites or in a pattern changing with every 2nd 2-hole supercell 
(Fig. [4]). The effect of isotopes on the thermal conductivity of 



graphene has recently attracted much attention due to the ad- 
vances in graphene synthesized by chemical vapour deposition 
(CVD). By using sequential input of 13 CH 4 and 12 CH 4 several 
groups have managed to fabricate graphene with regions of 
different isotopes lfT4ll . The local thermal conductivity of the 
regions with different isotopes has been measured by Raman 
spectroscopy lfT5l . The isotope effect has been addressed theo- 
retically for graphene fT6lL ifTTI and nanoribbons fT8ll . Inspired 
by this we here consider how the isotopes reduce the thermal 
conductance of GAL devices. In Fig. [4j\ we show T p h (scaled 
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Fig. 4. The effect of isotopes is considered for a long system (M = 12). 
The random isotopes of 13 C atoms are distribution at random positions at 
each M with a concentration of 50%. For the supercell pattern every second 
unit cell consist of 13 C isotopes. The temperature is 300K. 

by cell-width) for the longest system considered (M = 12 two- 
hole unit cells). We clearly see that isotope scattering give rise 
to a significant reduction in the thermal transmission. Compar- 
ing the reduction from the random isotope distribution to that 
of a pattern (alternating isotopes for every 2-hole supercell), 
it is clear that the "clustering" of isotopes leads to an extra 
reduction for especially the acoustic modes in the GAL region. 
Furthermore, we see that the difference will be less important 



at higher temperatures (> 300K used here). In Fig. [4j> we 
consider the thermal conductance for two GAL systems with 
different hole shape and size, as a function of length of GAL- 
device. For comparison, the thermal conductance is found to 
be 6.1nW/K for a pristine graphene device of the same width, 
and 0.41nW/K, 1.17nW/K, 0.76nW/K and 1.38nW/K for ideal 
{10,6arm}, {10,3arm}, {10,9zz} and {10,3zz} GALs, with 
length M > 5, respectively. With random isotopes the thermal 
conductance of the M = 1 2- {10,6arm} -device is reduced by 
55%, compared to 88% with the ordered isotope pattern. For 
devices of longer lengths we enter the diffusive region and 
cannot neglect phonon-phonon scattering as considered here. 

Since isotopes only affect the thermal properties we get 
a significant increase in ZT, especially at low gate bias 
where thermal transport is dominated by phonons. A peak 
max(ZT) « 0.8 is observed for the supercell pattern of 
isotopes at 300K, compared to max(ZT) « 0.26, 0.17, 0.07 
and 0.13 for the ideal {10,6arm}, {10,3arm}, {10,9zz} and 
{10,3zz} GALs without isotopes. The oscillations in ZT 
with M stems from the choice of patterning where we only 
have 50 — 50% isotopes for even M. With random isotopes 
we can exceed max(ZT) « 0.4 for the {10,6arm} lattice. 
This illustrates that one can benefit from a further reduction 
of the heat conductance in GAL devices to increase the 
thermoelectric efficiency. 

III. A MOLECULAR DYNAMICS APPROACH TO QUANTUM 
THERMAL TRANSPORT 

Finally, we will briefly discuss how the elastic quantum 
phonon transport can be obtained from generalized Langevin 
molecular dynamics(MD), as an alternative to the Green's 
function method. The method was proposed by Wang et 
al tTfj . and is a promising method for handling both low 
temperature quantum freeze-out of phonons, as well as the 
anharmonic interaction beyond perturbation theory. Both of 
these effects need to be included for graphene at 300K for 
devices longer than « 600A l28l . Anharmonic effects are 
clearly a very challenging problem for larger nanostructures 
such as GALs, where the role of the nanostructuring on the 
thermal properties needs clarification. 

We solve the generalized Langevin equation for the mass- 
normalized displacements, 
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by using the velocity-Verlet algorithm. The forces on the right 
hand side are from the full Brenner potential, the frictional 
damping force and the quantum noise. The two latter are 
due to the connection of the device to the two thermal 
baths (each lead). The essential variable in quantum thermal 
transport is the coupling matrix obtained from the retarded 
phonon self energy, Tb(u) = i(U r (uj) — II a (co>)). This object 
describes how much heat is absorbed(emitted) to(from) the 
surroundings. The non-Markovian damping kernel is obtained 
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Where r = £ — t' is the relative time so that the frictional 
damping force includes the memory at time t of a displacement 
at time t' . We choose a discrete time step At for a iVfc step long 
MD simulation (with corresponding Auj and discrete ujk in 
frequency space). To improve convergence and computation- 
time we include an artificial damping, e a d, to reduce the 
memory kernel length, n cut . We can from this calculate an 
effective imaginary part of the phonon self energy including 
this extra damping, 
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From this we can calculate the noise spectral power density 
from the quantum fluctuation-dissipation theorem: 
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We have implemented this method of generating a colored 
noise sequence and show several tests of the outcome with a 
large artificial damping below. The aim is to find out how 
accurate one has to represent the self-energy in terms of 
memory length and artificial damping to get results consistent 
with NEGF calculations. In Fig. [5] we show the imaginary 
part of the frequency dependent phonon bath self energy, and 
the noise correlation function. The first two figures (A and B) 
show how, even if we shorten the memory length (which in 
principle is as long as the simulation itself) to 100 time steps 
and increase the artificial damping to O.leV, the self energy 
is still approximated quite well. In the insert in Fig. |5p we 
have calculated the noise correlation function (noise product at 
two different times averaged over 10 4 starting points), which 
classically should equal kBT<f>(t — t'). This is not the case 
for graphene at 300K. As can be seen from the figure one 
has to go above 1000K for the classical fluctuation-dissipation 
theorem to be a good approximation. The reason for this is 
the quantum fluctuations, which give rise to a larger noise and 
kinetic energy below the Debye temperature, and the Bose- 
distribution of the phonons. Furthermore, we see from the last 
figure that MD with quantum heat baths gives results for the 
thermal conductance consistent with NEGF for short devices 
even at low temperatures. 

IV. Conclusion 

We have considered electronic and thermal transport prop- 
erties of finite GAL structures with disorder. We conclude that 
disorder can be a limiting factor for thermoelectric applications 
of GALs. A high degree of disorder either due to lattice 
irregularity or the lack of pure edge chirality will degrade 
the overall thermoelectric performance. However, the trend 
that an increased antidot density enhances the thermoelectric 
efficiency still survives. The presence of isotopes is found to 
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Fig. 5. A: Effect of memory length in the effective self energy used in 
MD. B: Effect of artificial damping in the effective self energy. C insert: 
Temperature dependence of the noise correlation function compared to the 
damping kernel. C: Thermal conductance from nonequalibrium MD with a 
temperature difference of 10% compared to the NEGF result for a {10,3arm} 
GAL with M = 1. e ad = 0.05eV, St = 0.5fs and n cut = 350 is used. 



drastically reduce the thermal conductance for GALs. Both 
in the case of lattice disorder and isotope mass- variation we 
have shown that disorder can be used to greatly enhance 
thermoelectric performance. In particular, if lattice disorder 
is kept at a minimum, and isotopes are changed in patterns on 
a length scale that affect most phonon modes, one can benefit 
from a very large reduction of the thermal conductance in such 
devices to obtain a good thermoelectric efficiency. Finally, we 
have illustrated how the harmonic quantum thermal transport 
can be obtained from generalized Langevin-MD simulations 
with quantum heat baths, yielding results similar to those 
obtained with the Green 's-function/Landauer approach. The 
MD approach is, in principle, able to address the effects of an- 
harmonicity in thermal transport of nano structured graphene. 
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